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SYSTEMATIC  ANALYSIS  OF  SATURN  WIRE  ARRAY 
IMPLOSION  OBSERVATIONS 


Introduction 

In  this  report,  we  will  describe  the  analysis  of  data  from  a  set  of  aluminum  wire-array  implosions 
conducted  on  the  Saturn  accelerator  at  Sandia  National  Laboratories.^’^  The  Saturn  experiments 
were  conducted  to  explore  the  validity  of  theoretical  implosion  models  by  investigating  the  variation 
of  the  keV  x-ray  yield  as  a  tradeoff  is  made  between  the  amount  of  mass  imploded  and  the  maximum 
velocity  achieved  in  the  implosion.  In  this  series  of  experimental  shots,  the  array  mass  per  unit 
length  m  and  the  initial  array  radius  tq  were  varied  in  a  way  designed  to  keep  the  total  kinetic 
energy  in  the  implosion  roughly  constant.  Thus,  variations  in  radiative  yield  should  reflect  only 
how  efficiently  this  kinetic  energy  is  converted  to  keV  radiation. 

In  the  experiments,  both  total  and  keV  radiative  yield  showed  a  discernible  trend  as  as  a  function 
of  either  array  mass  or  initial  radius.  One  goal  of  the  present  analysis  was  to  explain  this  trend 
by  gaining  greater  understanding  of  the  average  plasma  conditions  during  the  implosion.  Two 
important  factors  which  predict  output  power  in  the  present  experiments  are  the  average  ion  kinetic 
energy  at  stagnation  and  the  fraction  of  the  imploding  mass  involved  in  the  radiation. 

The  importance  to  K-shell  radiation  of  having  the  proper  ion  kinetic  energy  at  stagnation  has 
become  well  recognized.^  This  energy  must  be  large  enough  but  not  too  large:  The  implosion  must 
ionize  the  imploding  mass  and  heat  the  electrons  sufficiently  to  radiate,  but  should  not  cause  other 
processes  (stripping,  runaway  generation)  to  dominate.  This  theory  is  well  supported  in  the  present 
experiments. 

The  degree  of  mass  participation  in  the  radiation  also  strongly  influences  the  keV  output.^  In 
the  present  experiments,  we  found  considerable  variation  from  shot  to  shot  in  the  fraction  of  the 
imploding  mass  involved  in  K-shell  radiation  (from  0.04%  to  78%).  These  values  correlated  well 
with  the  K-shell  radiative  power. 

Experimental  observations  from  the  runs  consisted  of 

1 .  Filtered  PCD  output  records  as  a  function  of  time; 

2.  Time-integrated  keV  spectra; 

3.  Pinhole  photographs,  both  time-integrated  and  at  time  intervals. 

The  systematic  analysis  given  here  uses  a  procedure  developed  earlier  at  NRL.'*  Using  this 
procedure,  we  can  infer  from  x-ray  photodiode,  spectroscopic  and  pinhole  measurements  the 
average  radiated  power,  the  average  ion  density  and  plasma  temperature,  and  the  average  fraction 
of  total  mass  that  was  involved  in  K  shell  x-ray  radiation. 

This  report  is  organized  as  follows.  First,  a  simple  implosion  model  is  presented,  to  introduce 
the  terminology  and  basic  scaling  behavior.  Parameters  for  the  experimental  runs  are  presented. 
Next,  the  details  are  given  for  analysis  of  the  pulsewidth,  line  spectrum,  and  pinhole  photograph 
data.  The  measurements  from  the  data  are  then  summarized  and  trends  are  pointed  out. 
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Simple  implosion  model 


The  slug  model  gives  a  good  basic  framework  for  discussion  of  the  implosion  dynamics.  In 
this  model,  all  of  the  imploding  material  forms  a  cylindrically  symmetric  shell  at  a  radius  r,  which 
obeys  the  force  equation 

(fr  \dL 

7Tl~- —  — - J  /i\ 

dt^  21  dr 

where  m  is  the  mass  per  unit  length,  I  is  the  wire  length,  I{t)  is  the  total  current,  and  L{t)  is  the 

inductance  of  the  circuit  element  consisting  of  the  imploding  material  and  the  concentric,  external 
anode.  For  an  array  of  N  wires, 

^  ^  N-  \  I  ^ 
dr  N  27r  r  r 

where  /Xq  =  47r  x  10~^  H  m~^. 

The  current  through  the  pinch  at  a  given  time  is  determined  by  the  parameters  of  the  transmission 
line  circuit,  which  has  a  total  inductance  Lt  and  impedance  Zo.  The  circuit  equation  is; 

dLrl  , 

—  +IZo  =  V(t)  (3) 

The  simultaneous  solution  of  Eq.  (1)  and  Eq.  (3)  for  a  given  voltage  waveform  F(f)  gives  the  radial 
position  and  velocity  of  the  implosion,  as  long  as  plasma  pressure  is  negligible. 

In  designing  the  parameters  for  the  experiments  described  here,  both  the  force  and  circuit 
equations  were  solved  for  an  input  V(t)  appropriate  for  Saturn.  Useful  insight,  however,  can  also 

be  gained  by  neglecting  the  circuit  equation  and  assuming  a  simplified  linear  time  dependence  for 
the  current:^ 

I(t)  =  lit.  (-4) 

To  aid  in  the  description  of  the  experimental  results,  we  will  summarize  results  obtained  earlier^ 
for  the  solution  of  Eq.  (1)  with  I(t)  from  Eq.  (4). 

First,  we  write  Eq.  (1)  in  dimensionless  form,  by  defining 

x  =  r/ro;  T  =  t!to;  i  =  I/Io  (5) 

where  ro  is  the  initial  array  radius,  and  the  constants  to  and  lo  are  left  arbitrary  for  now.  The 
dimensionless  form  of  Eq.  (1)  is 


(6) 

miroltoY  ■ 

(7) 

2 


Because  of  Eq.  (4),  current  and  time  are  in  one-to-one  correspondence.  We  can  combine  these  by 
defining 


q=  [A^^\hto/Iof'^]  T 


(8) 


so  that  Eq.  (6)  becomes 


(Px  , 
X- —  =  —q 

dq^ 


(9) 


with  the  boundary  conditions  1(0)  =  1  and  q(0)  =  0. 

Eq.  (9)  can  be  numerically  solved  to  give  a  curve  of  x(q),  which  is  the  general  radial  trajectory 
of  the  slug-model  implosion  with  a  linear  current  rise.  This  solution  and  that  for  dx  jdq  are  shown 
in  Fig.  (1).  The  actual  radial  values  are  given  by  r  =  Tox{q{t)),  with  the  appropriate  conversion 
factors: 


q(t)  = 


mr^ 


(10) 


Note  that  this  relation  is  independent  of  the  choice  of  the  normalization  constants  Iq  and  to,  which 
have  not  yet  been  defined.  Characteristic  values  that  might  be  used  for  normalization  are  the 
maximum  current  Ijnax  =  htmax  and  the  circuit  relaxation  time  Ilr  =  Lo/Zq. 

In  this  model,  all  implosions  behave  in  the  same  way,  following  the  curve  in  Fig.  (1)  from 
a:(0)  =  1  to  stagnation  at  x{qf)  =  Xf.  Stagnation  occurs  at  the  point  where  the  plasma  pressure, 
omitted  from  Eq.  ( 1 ),  exceeds  the  inward  magnetic  pressure.  We  assume  that  this  pressure  becomes 
important  too  suddenly  to  affect  the  overall  dynamics  (i.e.,  the  zero-pressure  curve  x{q)),  and 
arbitrarily  define  the  stagnation  position  Xf  based  on  experimental  observations.  The  usual 
definition  has  been^  for  0.1  <  Xf  <  0.2;  this  choice  is  supported  by  the  pinhole  photograph 
data  presented  here.  The  stagnation  time  tf  is  derived  from  g/  by  using  Eq.  (10).  The  stagnation 
point  is  also  the  time  of  maximum  implosion  velocity  (proportional  to  dxjdq). 

An  important  implosion  parameter  is  the  kinetic  energy  per  ion^  Ki  =  ^rriivj,  where  m;  is  the 
ion  mass  and  Vf  is  the  ion  velocity  at  the  stagnation  time  tf.  This  can  be  written  as 

T  2 


Ki  =  fii 


(  dx\ 


(11) 


where  the  energy  conversion  factor  fii  =  1.398x  10^  eV  cm~^  ns.  For  efficient  K-shell  emission,  Ki 
must  exceeed  an  energy  Emin,  which  is  required  to  strip  all  but  one  or  two  electrons  from  the  atom 
and  heat  the  stripped  electrons  sufficiently  to  radiate.  For  aluminum,  Emin  is  roughly  12  keV/ion.^’^ 
The  parameter 

ri  =  KilEmin  (12) 


indicates  the  degree  to  which  the  minimum-energy  condition  is  fulfilled.  We  expect  that  K-shell 
radiation  will  be  greatest  when  q  >  I  but  not  forq  >>  1,  when  processes  other  than  K-shell 
radiation  will  dominate. 
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The  energy  parameter  determining  the  maximum  total  radiation  is  the  kinetic  energy  per  unit 
length^  Kc  =  which  is  given  by 


If 


(13) 


We  expect  that  the  total  x-ray  yield  will  increase  as  is  increased  (though  was  held  constant 
in  the  experiments  analyzed  here). 

The  machine  parameters  affecting  the  implosion  are  mass  m,  initial  radius  tq  and  current  rise 
rate  7i  (or  maximum  current  I^ax  =  htf)-  Since  qf  and  (dx/dq)f  are  the  same  for  all  implosions, 
the  following  scaling  must  apply  for  the  implosion  time,  particle  energy  (for  a  given  element),  and 
total  energy: 


tf  ~  [mrl]^/yirnax  (14a) 

V'-rn-^lLx  (146) 

Kc  ~  I^ax-  (14c) 

The  77  value  is  independent  of  tq  and  Kc  is  independent  of  both  m  and  tq. 

The  shot  parameters  for  the  two  experimental  series  are  given  in  Table  I.  Two  sets  of  6  shots 
were  made,  one  set  using  larger  initial  radii  (and  anode  diameter)  than  the  other.  The  implosion 
time  for  these  shots  varies  from  about  50  ns  at  low  mass  to  90  ns  at  high  mass,  based  on  the  linearly 
rising  current  model;  with  the  Saturn  circuit  model,  these  times  would  be  slightly  shorter. 

The  accuracy  of  the  simple  model  in  determining  implosion  energy  can  be  assessed  by 
comparing  the  tabulated  values  of  the  energy  parameter  q,  obtained  using  the  more  accurate 
Saturn  circuit  model  (Eq.  (3)),  and  of  770.2,  the  same  parameter  computed  from  the  linear  current 
model  with  Xf  =  0.2.  The  difference  is  generally  under  10%.  The  experimental  design  of  holding 
Kc  constant  is  reflected  in  the  invariance  of  the  product  77777  =  2x  10^  The  successful  achievement 
of  this  goal  can  be  seen  in  the  near  constancy  of  the  maximum  current 

Pulsewidth  analysis 

The  pulsewidth  was  evaluated  from  the  net  keV  and  total  radiation  observations.  Output  traces 
from  the  filtered  x-ray  PCD  s,  which  had  a  resolution  of  better  than  10  ns,  are  shown  in  Fig.  (2). 
The  pulsewidth  was  defined  to  be  the  full  width  at  half-maximum,  which  quantity  was  uniquely 
given  in  the  present  data.  Note,  however,  that  two  pulses  were  often  present,  a  strong  initial  pulse 
and  decay  followed  by  a  “bounce”  pulse.  The  bounce  is  prominent  in  the  high-mass  shot  1498  and 
is  present  in  all  of  the  large  initial  radius  shots  (set  2,  Fig.  (2b)).  In  these  observations,  however, 
the  half-maximum  width  never  extended  beyond  the  first  (and  stronger)  peak. 
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Table  I.  Shot  parameters 


Shot 

Array  dia 

Mass 

^max 

Total  yield 

keV  yield 

Vo.2 

V 

No. 

(cm) 

ifig/cm) 

(MA) 

(kJ) 

(kJ) 

Set  1:  Anode  diameter  2.80  cm 

1483 

1.71 

328 

8.9 

224 

76.0 

6.9 

6.0 

1484 

1.59 

473 

8.7 

242 

38.1 

5.7 

4.2 

1487 

1.46 

739 

9.0 

291 

30.0 

3.1 

2.7 

1488 

1.40 

949 

9.4 

224 

20.0 

2.7 

2.2 

1489 

1.34 

1313 

8.8 

270 

13.2 

1.8 

1.6 

1498 

1.28 

2052 

9.8 

228 

4.9 

1.3 

0.95 

Set  2:  Anode  diameter  4.80  cm 

1492 

3.04 

118 

9.1 

172 

8.0 

20 

22 

1493 

2.81 

161 

9.6 

170 

12.2 

16 

17 

1494 

2.73 

185 

9.4 

196 

22.4 

14 

15 

1495 

2.42 

328 

9.7 

331 

63.3 

8.2 

8.3 

1496 

2.29 

473 

9.7 

364 

61.2 

5.7 

5.6 

1497 

2.20 

739 

9.9 

338 

30.1 

3.8 

3.5 

Results  for  the  keV  pulsewidth  are  plotted  in  Fig.  (3),  as  a  function  of  both  array  radius  and 
mass.  While  the  data  show  some  variability,  generally  the  pulsewidth  increases  with  mass  and 
decreases  with  increasing  initial  radius.  The  qualitative  differences  in  the  output  traces  of  Fig.  (2a) 
and  Fig.  (2b)  will  also  be  reflected  in  the  pinhole  data  presented  below,  with  broader,  more  irregular 
output  traces  associated  with  diffuser,  lower-contrast  pinhole  images. 

Spectral  analysis 

The  time-integrated  x-ray  spectra  of  the  two  sets  of  Saturn  shots  are  shown  in  Fig.  (4).  The 
important  K-shell  lines  are  listed  in  Table  II.  These  lines  can  be  clearly  seen  in  the  spectra  of 
almost  all  the  experimental  shots.  While  spectral  broadening  and  properties  of  satellite  lines  could 
in  principle  be  used  to  obtain  information  on  plasma  conditions,  experimental  conditions  dictate 
the  line  energy  as  the  most  reliable  spectral  diagnostic.  Peak  widths  depend  on  the  time-varying 
and  parameter-dependent  width  of  the  source,  and  so  provide  no  readily  obtainable  information 
on  plasma  parameters.  The  actual  power  emitted  in  a  single  line,  given  by  the  integral  under  the 
spectral  curve,  cannot  be  determined  with  certainty  because  of  the  difficulty  of  calibrating  the 
instrumentation,  and  of  accurately  determining  continuum  radiation.  The  ratio  between  power 
emitted  in  different  lines,  however,  is  independent  of  instrumental  calibration,  and  is  a  sensitive 
function  of  plasma  density  and  temperature.  We  use  the  ratio  of  the  most  prominent  lines,  the  Ly  a 
and  He  a  lines,  in  our  analysis. 

The  energy  radiated  in  each  line  was  found  by  computer  analyzing  the  digitized  spectra.  The 
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Table  II.  Energies  of  important  A1  K-shell  lines  and  edges 


Line 

Energy  (eV) 

Wavelength  (A) 

He  IC 

1588 

7.81 

He  a 

1598 

7.76 

Ly  a 

1728 

7.17 

He/3 

1869 

6.64 

He  7 

1964 

6.31 

He  6 

2006 

6.18 

He  edge 

2086 

5.94 

Ly/3 

2048 

6.05 

Ly7 

2160 

5.74 

Ly  8 

2212 

5.61 

H  edge 

2304 

5.38 

strategy  used  was  as  follows.  To  minimize  interference  from  data  jitter,  we  analyzed  not  the 
spectrum  itself,  but  its  integral.  These  integrals  are  shown  in  Fig.  (5),  and  are  sized  so  they  may 
be  overlaid  onto  the  corresponding  spectra  of  Fig.  (4).  Clearly,  a  spectral  line  causes  a  sudden 
increase  in  the  slope  of  the  integral  curve.  The  energy  radiated  in  the  line  is  the  increase  along  the 
ordinate  of  the  integrated  curve  from  where  its  slope  sharply  increases  to  where  its  slope  levels  off. 
The  magnitude  of  the  slope  of  the  integrated  curve  outside  the  the  spectral  lines  is  proportional  to 
the  baseline  height  of  the  spectrum,  or  to  the  amount  of  continuum  radiation  at  a  given  energy.  For 
the  present  analysis,  no  attempt  was  made  to  subtract  such  “background”  radiation  from  the  total 
line  energy. 

To  identify  spectrum  peaks,  the  average  or  baseline  slope  of  the  integral  curve  was  computed 
as  the  slope  between  the  initial  point  and  the  current  point.  This  was  compared  to  the  "current" 
slope  (average  of  slope  at  next  2  pts).  A  peak  was  deemed  to  begin  whenever 

(current  slope)  >  cq  *  (baseline  slope),  (15) 

where  Cq  «  2  is  a  convenient  factor,  whose  exact  value  is  not  important  for  the  sharp  peaks  found 
here.  Two  strategies  were  tried  for  numerically  defining  the  end  of  the  peak: 

1.  (current  slope)  <  Cj  *  (baseline  slope)  (16o) 

2.  (current  slope)  <  (max  slope  during  peak)/ Cl  (166) 

where  ei  «  2  is  another  arbitrary  factor  chosen  for  numerical  convenience.  The  criterion  in  Eq.  ( 1 6a) 
combines  sub-peaks,  and  is  insensitive  to  forked  peaks  or  satellites.  Thus,  if  the  continuum  signal 
increases  on  the  high-energy  side  line  peak  (typically  so  at  the  high-energy  end  of  spectrum),  this 
criterion  tends  to  miss  the  end  of  that  peak  and  so  group  all  peaks  in  the  tail  together.  The  criterion 
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in  Eq.  (16b)  sensitively  resolves  almost  all  peaks,  but  may  be  too  sensitive  when  there  is  jitter  in 
the  data.  The  first  strategy  was  used  for  the  present  results. 

Some  results  from  the  spectroscopic  analysis  are  shown  in  Table  III.  The  range  of  photon 
energies  between  1.5  keV  and  2.5  keV  was  used  in  the  analysis;  this  range  was  of  primary  interest 
and  was  also  where  the  instruments  were  most  accurate.  This  energy  range  typically  included 
50-60%  of  the  total  energy  in  the  obtained  spectra;  the  fraction  for  each  shot  is  given  in  the  second 
column  of  the  table.  The  third  column  shows  the  fraction  of  the  total  energy  in  the  restricted  range 
that  was  emitted  in  line  radiation.  Later  columns  show  the  fraction  of  the  total  restricted  energy  that 
was  emitted  in  each  of  the  lines  whose  energies  are  listed  in  Table  II;  these  should  sum  roughly  to 
the  total  in  the  third  column.  Line  satellites  are  generally  included  with  the  line  totals.  Line  ratios 
are  found  by  simply  taking  the  ratio  of  the  column  figures  for  the  desired  lines;  for  the  present 
analysis,  we  used  the  ratio  of  Ly  a  to  He  a  (+  intercombination)  lines,  since  these  were  always  the 
most  prominent. 


Table  HI.  Fraction  of  keV  radiated  energy  in  atomic  lines 


Shot 

No. 

Fraction  in 

1. 5-2.5  keV 

Lines 

He  a 

Fraction  of  restricted  spectrum  energy 

Ly  a  He  /?  He  7  He  ^  Ly  /? 

Ly  7 

Ly  6 

1483 

0.618 

0.164 

0.242 

0.035 

0.017 

0.017 

0.092 

0.024 

0.015 

1484 

mSSM 

0.593 

0.202 

0.253 

0.034 

0.017 

0.017 

0.028 

0.013 

0.009 

1487 

0.614 

0.199 

0.276 

0.034 

0.016 

0.018 

0.028 

0.015 

0.009 

1488 

mBM 

0.647 

0.226 

0.298 

0.031 

0.016 

0.013 

0.028 

0.013 

0.006 

1489 

0.664 

0.666 

0.241 

0.306 

0.033 

0.015 

0.014 

0.026 

0.012 

0.005 

1492 

0.106 

0.709 

0.291 

0.357 

0.056 

1493 

1.000 

0.560 

0.188 

0.238 

0.034 

0.011 

0.006 

0.053 

0.018 

0.007 

1494 

0.644 

0.821 

0.238 

0.375 

0.050 

0.016 

0.017 

0.079 

0.028 

0.011 

1495 

0.006 

0.522 

0.187 

0.256 

0.038 

0.017 

0.018 

1496 

0.555 

0.547 

0.200 

0.211 

0.027 

0.016 

0.016 

0.031 

0.016 

0.009 

1497 

0.517 

0.543 

0.232 

0.232 

0.035 

0.017 

0.015 

1498 

0.603 

0.712 

0.269 

0.338 

0.028 

0.012 

0.008 

0.023 

0.010 

0.003 

Pinhole  photograph  analysis 

The  main  goal  in  analysis  of  the  time-integrated  pinhole  data  was  to  determine  the  width  of 
the  keV  radiating  region.  Two  time-integrated  photographs  were  taken,  one  high-gain  and  one 
low-gain,  as  well  as  10  fixed-time  snapshots,  taken  at  intervals  of  3  ns  for  a  1  ns  exposure  time. 
Data  acquisition  was  electronic,  and  corrections  were  made  for  the  detector  response  before  data 
analysis  was  begun. 
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The  high-gain,  time-integrated  photograph  was  used  for  the  present  analysis.  This  image 
was  entered  into  the  computer,  and  an  algorithm  was  developed  to  obtain  an  average  value  of  the 
pinch  width  that  correlated  well  with  visual  observation,  even  for  irregular  implosions  and  faint 
images.  No  densitometry-type  calibration  was  made,  because  the  degree  of  arbitrariness  inherent 
in  the  determination  of  pinch  radius  for  any  one  shot  made  such  precision  unnecessary.  However, 
because  the  same  determination  procedure  was  consistently  followed  for  all  shots  in  a  given  run, 
the  pinhole  analysis  will  accurately  identify  the  trends  in  the  data. 

When  received  for  analysis,  the  pinhole  images  were  on  film.  They  were  then  digitized,  using 
a  Hewlett-Packard  ScanJet  II  scanner  and  accompanying  software.  Software  settings  were  150 
dpi,  256  black  and  white  shades,  brightness  120  and  contrast  131.  The  choice  of  resolution  and 
color  bit  length  was  made  to  optimize  information  content  while  keeping  the  file  size  reasonable. 
A  full-page  scan  of  all  12  images  for  a  given  shot  typically  generated  a  TIFF  file  of  about  500K 
bytes.  One  pinhole  image  extracted  from  this  file  takes  up  about  41K  bytes. 

One  possible  algorithm  for  finding  the  implosion  width  would  be  to  average  the  pixel  values 
in  the  z  direction,  taking  the  width  of  the  radiating  region  to  be  the  width  at  the  half-maximum 
or  similar  pixel  value.  The  most  important  quantity  for  width  determination  is  the  variation  in 
intensity  from  the  peak,  not  the  absolute  intensity,  which  is  more  a  function  of  the  overall  power 
of  the  implosion.  An  illustration  of  the  method  of  averaging  over  z  is  shown  in  Fig.  (6),  where 
pinhole  photographs  are  shown  with  the  r-dependent  averaged  pixel  values.  For  both  cases  shown, 
the  z-averaged  film  darkens  progressively  until  a  point  of  maximum  darkness  is  reached.  This 
maximum  point  is  taken  to  be  the  pinch  center.  The  pinch  width  can  be  obtained  by  defining  some 
fraction  of  the  peak  value  as  a  criterion  (e.g.,  half-maximum).  This  simple  procedure  works  best 
with  well-defined,  straight  implosions  like  that  in  shot  1484.  Lower-contrast,  kinked  implosions 
like  that  of  shot  1498  are  more  difficult  to  analyze. 

The  problem  with  the  z-averaging  method  just  described  is  that  when  the  implosions  are  not 
straight,  or  when  the  contrast  is  so  low  that  there  are  multiple  peaks  in  the  averaged  profile, 
the  averaging  method  gives  unrealistically  large  pinch  widths.  To  follow  a  crooked  or  kinked 
implosion,  it  is  better  to  conduct  the  2  average  over  multiple  small  segments,  rather  than  over  the 
whole  pinch.  After  some  trial,  it  was  found  that  a  6-pixel  long  segment  was  optimum:  for  shorter 
segment  lengths,  the  radial  profile  became  erratic,  while  for  larger  segment  lengths,  the  computed 
widths  became  larger,  but  at  a  very  slow  rate. 

The  results  of  this  analysis  for  the  two  sets  of  implosions  are  shown  in  Fig.  (7)  and  Fig.  (8).  In 
these  figures  are  shown,  as  a  function  of  z,  the  computed  width,  the  contrast  in  the  image  (possible 
pixel  values  run  from  0  to  255),  and  the  position  of  the  center  point  (the  darkest  point  in  the  image). 
To  obtain  the  pinch  width,  the  average  was  taken  of  the  widths  at  each  point  in  2.  These  widths  were 
determined  just  as  was  shown  in  Fig.  (6),  except  that  instead  of  averaging  over  the  whole  pinch 
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length,  a  6-pixel  average  was  taken.  The  total  width  included  only  those  points  with  sufficiently 
high  image  contrast  that  the  computed  width  was  reliable  (contrast  >  50);  in  low-contrast  areas, 
the  width  and  center  position  are  primarily  a  result  of  random  noise. 

The  procedure  used  to  determine  the  width  of  the  radiating  region,  then,  was  to  scan  the 
pinhole  images,  extract  the  high-gain  time-integrated  pinhole  image,  average  the  pixel  values  over 
all  6-pixel  long  axial  segments  in  the  image,  and  define  the  width  as  the  average  of  the  full-width 
quarter-maximum  (FWQM)  values  at  each  point  in  z  where  there  was  sufficient  contrast  to  discern 
an  image.  The  method  gave  consistent  and  reasonable  values  for  the  width  for  almost  all  of  the 
pinhole  images,  even  when  the  image  contrast  was  too  low  for  reliable  evaluation  by  eye.  Shots 
where  there  was  a  problem  determining  width  were  1493,  which  was  very  low  contrast;  1496, 
where  evident  light  contamination  ruined  the  shot  image;  and  1487,  for  which  no  pinhole  images 
where  available.  In  addition,  no  pinhole  image  was  available  for  shot  1483,  but  the  pinhole  image 
from  another  shot  with  similar  parameters  and  yield,  shot  1606,  was  substituted. 

The  size  for  each  shot  of  the  radiating  region,  derived  from  the  pinhole  photograph 
measurements,  is  plotted  in  Fig.  (9),  as  a  function  of  both  mass  and  array  diameter.  Generally, 
the  radiating  region  contracts  as  the  mass  increases  and  the  initial  array  diameter  decreases.  The 
exceptions  to  this  rule  are  the  two  low-mass,  large-diameter  shots  (1492  and  1493),  where  in  any 
event  the  pinhole  images  were  so  faint  as  to  make  the  diameter  measurements  suspect. 

In  principle,  the  length  of  the  radiating  region  could  be  determined  with  the  same  procedure  as 
was  used  for  the  width.  In  the  pinhole  photographs  used,  however,  obstructions  in  the  field  of  view 
made  this  impossible,  and  it  was  assumed  that  in  every  case  there  was  K-shell  radiation  over  the 
full  2  cm  length  of  the  pinch.  The  photographs  show  that  this  is  a  reasonable  assumption  at  least 
in  the  field  of  view,  for  hot-spot  formation,  while  occasionally  evident,  did  not  significantly  reduce 
the  axial  extent  of  the  averaged  radiating  region. 

Results 

Data  interpretation  followed  the  procedure  of  Coulter,  Whitney  and  Thornhill.^  Basically,  this 
procedure  is  to  use  the  measurements  just  described  and  compare  them  with  theoretical  predictions 
for  a  range  of  densities  and  temperatures.  Contours  can  be  drawn  in  density  and  temperature  for 
the  observed  power  per  unit  length  as  well  as  for  the  observed  line  ratios.  The  place  where  these 
contours  intersect  marks  the  inferred  average  plasma  density  and  temperature. 

The  theoretical  calculations  were  performed  for  a  uniformly  dense,  isothermal  plasma  of  the 
measured  diameter  and  2  cm  in  length,  which  was  in  collisional-radiative  equilibrium  (CRE).® 
Variation  in  opacity  with  size  of  the  emission  region  was  included  in  the  calculation.  Theoretical 
contours  of  keV  power  emitted  per  unit  length  are  shown  in  Fig.  (10),  for  a  plasma  diameter  of 
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3  mm.  The  experimental  power  per  unit  length  was  obtained  using  the  yield  from  Table  I  and  the 
measured  pulsewidth;  this  selects  one  of  the  contours  in  Fig.  (10).  The  ratio  of  Ly  a  to  He  a  (plus 
intercombination  line)  power  was  calculated  in  CRE/  and  contours  are  shown  in  Fig.  (11).  The 
experimental  line  ratio,  obtained  from  Table  III,  selects  the  appropriate  contour  in  Fig.  (11).  The 
intersection  of  the  two  selected  contours  marks  the  average  plasma  and  density. 

For  example,  in  shot  1483  the  total  extent  of  the  radiating  region  is  3.1  mm,  the  average  keV 
power  is  4.2  TW  (2.1  TW/cm)  and  the  line  ratio  is  1.48.  The  intersection  of  the  marked  contours 
is  at  an  electron  temperature  of  600  eV  and  an  ion  density  of  7  x  10^®  cm~^. 

Once  the  average  density  has  been  determined,  the  fraction  of  the  mass  involved  in  keV 
radiation  can  be  estimated  by  comparing  the  mass  in  the  radiating  region  (whose  extent  has  been 
measured)  with  the  total  mass  in  the  wire  array.  The  fraction  obtained  in  this  manner  might  be 
inaccurate  because  of  deviations  from  cylindrical  symmetry  (tending  to  make  the  calculation  an 
overestimate)  or  because  of  radial  density  variations  (tending  to  make  it  an  underestimate).  But 
this  fraction  does  give  some  idea  of  the  effectiveness  of  the  machine,  for  the  given  parameters,  in 
including  all  available  mass  in  the  radiation  pulse. 

The  measurements  from  the  data,  and  the  plasma  parameters  they  imply,  are  summarized  in 
Table  IV.  The  widths  for  both  the  total  radiation  pulse  and  for  the  keV  radiation  pulse  Tk,v 
are  given,  with  the  corresponding  average  power  given  by  the  energy  yield  (from  Table  I)  divided 
by  the  pulsewidth.  The  line  ratio  is  taken  from  Table  III,  with  (as  in  that  table)  the  He  a  fraction 
including  the  intercombination  line.  The  average  diameter  of  the  radiating  region  is  taken  to  be 
the  FWQM  of  the  6-pixel  z-averaged  digitized  pinhole  image,  as  described  above.  The  effective 
energy  parameter,  for  the  implosion  is  an  estimate,  based  on  the  linear  current  model  described 

above  and  observed  stagnation  position  Xf  =  (FWQM  radius) /ro.  Those  parameters  are  used  in 
Eq.  (11)  to  get  the  ion  kinetic  energy  and  t]  value.  The  average  ion  density,  electron  temperature, 
and  fraction  of  the  mass  radiating  in  the  K  shell  are  obtained  from  the  data  measurements  as  just 
described. 

From  Table  IV,  it  is  evident  that  both  the  total  and  keV  radiated  power  show  a  regular  variation 
with  mass  and  array  diameter.  The  keV  power  is  plotted  against  these  parameters  in  Fig.  (12).  The 
variations  in  power  are  not  explained  by  variations  in  the  average  density  and  temperature,  which 
are  shown  in  Fig.  (13)  and  Fig.  (14).  A  similar  pattern  is  found,  however,  in  the  radiating  mass 
fraction,  shown  in  Fig.  (15).  The  participating  mass  fraction  was  evidently  a  very  important  factor 
in  determining  the  keV  output  power  in  these  runs.  The  keV  power  is  plotted  as  a  function  of  mass 
fraction  in  Fig.  (16);  although  the  relationship  is  not  monotonic,  a  greater  mass  fraction  generally 
results  in  greater  power  output.  The  effect  on  the  keV  power  of  varying  rj  is  also  shown  in  Fig.  (16). 
As  was  expected,  the  greatest  keV  output  came  from  the  moderate-77  shots,  with  maximum  output 
for  6  <  7/  <  8.  The  same  trend  is  seen  with  the  calculated  value  T^e//,  except  that  the  keV  power 
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Table  FV.  Results  of  data  analysis 


Shot 

Total  power 

TkeV 

keV  power 

Lya/Hea 

FWQM  dia 

Vcff 

Tli 

Te 

K-shell 

No. 

(ns) 

(TW) 

(ns) 

(TW) 

ratio 

(cm) 

(10 

19  cm-3) 

(eV) 

fraction 

1483 

20 

11.20 

18 

4.22 

1.48 

0.31 

13 

7.0 

600 

0.72 

1484 

25 

9.68 

22 

1.73 

1.26 

0.21 

lA 

7.0 

540 

0.23 

1487 

25 

11.64 

20 

1.50 

1.39 

0.20(“) 

3.9(“) 

6.1 

600 

0.12 

1488 

30 

7.47 

23 

0.87 

1.32 

0.19 

3.2 

4.5 

620 

0.06 

1489 

51 

5.29 

31 

0.43 

1.26 

0.19 

2.8 

2.0 

830 

0.02 

1498 

87 

2.62 

32 

0.15 

1.25 

0.09 

2.7 

2.7 

760 

0.004 

1492 

29 

5.93 

25 

0.32 

1.24 

0.12 

43.2 

3.5 

700 

0.15 

1493 

26 

6.54 

27 

0.45 

1.27 

0.24 

25.4 

1.7 

900 

0.21 

1494 

27 

7.26 

26 

0.86 

1.58 

0.45 

15.7 

1.2 

1150 

0.46 

1495 

22 

15.05 

19 

3.33 

1.37 

0.41 

9.1 

4.3 

620 

0.78 

1496 

36 

10.11 

32 

1.91 

1.06 

0.32('’) 

7.1 

5.0 

500 

0.38 

1497 

38 

8.89 

36 

0.84 

1.00 

0.28 

4.9 

2.8 

600 

0.10 

No  pinhole  data;  value  inferred  from  trend 
Unreliable  data:  Actual  image  is  washed  out 


maximum  occurs  at  T/e/y  «  10. 

Conclusions 

In  these  experiments,  the  total  implosion  kinetic  energy  was  held  constant,  while  mass  and  initial 
radius  were  varied  over  a  wide  range.  A  detailed  and  consistent  analysis  of  the  experimental  data 
was  conducted  to  estimate  the  average  plasma  parameters  during  the  implosion,  and  to  investigate 
trends  in  these  that  could  explain  the  variation  in  keV  radiated  power. 

As  is  shown  in  Fig.  (16),  both  the  average  energy  per  particle,  represented  by  rj,  and  the 
fraction  of  the  array  mass  participating  in  the  implosion  showed  a  clear  relationship  to  the  keV 
radiated  power.  No  such  simple  relationship  was  found  with  other  plasma  parameters.  The 
consideration  of  mass  participation  could  be  particularly  important  on  high-power,  short-pulsewidth 
machines  like  Saturn.^  With  a  short  pulsewidth,  mass  participation  could  be  reduced  by  incomplete 
wire  vaporization,  and  initial  nonuniformities  in  plasma  temperature  and  density  will  be  more 
pronounced,  affecting  current  penetration  and  implosion  dynamics.  At  high  power  levels,  both 
macro-  and  micro-instabilities  may  complicate  the  implosion  dynamics,  and  L-  and  M-shell  losses 
could  become  significant  during  the  implosion. 

The  results  presented  here  can  only  provide  information  on  the  “average”  state  of  the  plasma 
during  the  implosion.  Lacking  data  with  significant  time  and  spatial  resolution,  we  could  not 
consider  the  significant  time  variations  in  plasma  conditions  and  radiation;  nor  could  we  include  the 
significant  gradients  in  temperature,  density,  and  electric  and  magnetic  fields  that  occur  during  the 
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irnplosion.  Still,  these  self-consistent  calculations  provide  valuable  insight  into  the  average  plasma 
conditions  during  the  radiation  pulse,  and  also  provide  a  consistent  means  of  comparing  different 
PRS  implosions.  To  understand  the  impact  of  dynamic  variations  and  spatial  gradients,  MHD 
simulations^  can  be  used  to  generate  data  for  known  implosion  parameters  that  can  be  analyzed 
with  the  same  procedures  given  here  (such  an  effort  is  currently  underway).  As  temporally-  and 
spatially-resolved  experimental  measurements  become  available,  the  procedures  detailed  here  can 
be  readily  used  to  provide  rapid  analysis  of  the  data. 
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X  (normalized  radius) 


Implosion  position  and  velocity 
with  Linear  current  rise 


-  position 

- velocity 
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q  (normalized  time) 


Figure  1.  Solution  to  Eq.  (9)  for  normalized  radial  position  x  =  r/ro  and  normalized  implosion 
velocity  dxfdq,  as  a  fimction  of  normalized  time  9,  given  by  Eq.  (10).  In  the  linear-current  model, 
every  implosion  follows  these  curves  until  x  =  x/,  where  x/  is  defined  to  be  the  stagnation  point. 
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dx/dq  (normalized  implosion  velocity) 


PCD  OUTPUT  -  SET  1 


•tr  180 


(0  120 

c 

D3 


190  230  270  310  350 


Time  (ns) 


190  230  270 


Time  (ns) 


XJ  120 

k— 

ca 


150  190  230 


Time  (ns) 


Time  (ns) 


150  190  230  27 

Time  (ns) 


Time  (ns) 


Figure  2a.  X-ray  photodetector  output  for  each  experimental  shot,  used  to  determine  the  x-ray  pulsewidth 
The  total  signal  is  given  by  the  dashed  lines,  and  the  filtered  keV  signal  is  given  by  solid  lines. 

The  set  grouping  is  given  in  Table  I. 


keV  pulsewidth 


Mass  loading  (^ig/cm) 
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Figxire  3.  The  keV  pulsewidth,  measured  from  Fig.  (2)  as  described  in  the  text,  plotted  as  a  function 
of  array  linear  mass  density  and  initial  array  diameter.  Points  from  data  set  1  (see  Table  I)  are 
plotted  with  a  plus  sign  (+),  and  points  from  data  set  2  are  plotted  with  a  star  (*). 
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Figure  4a.  Measured  spectra  for  each  experimental  shot.  Important  lines  (listed  in  Table  H)  are  identified. 
The  set  grouping  is  given  in  Table  I. 
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Figure  4b.  Measured  spectra  for  each  experimental  shot.  Important  lines  (listed  in  Table  H)  are  identified. 
The  set  grouping  is  given  in  Table  I. 


18 


SET  1  INTEGRATED  SPECTRA 
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Figure  5a.  Integral  under  the  curve  of  each  of  the  spectra  in  Fig.  (4),  and  used  in  the  data  analysis.  The  power 
emitted  in  each  spectral  line  is  proportional  to  the  rise  in  the  integral  curve  over  the  width  of  that  line.  Prominent 
peaks  in  Fig.  (4)  are  associated  with  a  steep  slope  in  the  integral  curves.  The  set  grouping  is  given  in  Table  I. 
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SET  2  INTEGRATED  SPECTRA 
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Figure  5b.  Integral  under  the  curve  of  each  of  the  spectra  in  Fig.  (4b),  and  used  in  the  data  analysis.  The  power 
emitted  in  each  spectral  line  is  proportional  to  the  rise  in  the  integral  curve  over  the  width  of  that  line.  Prominent 
peaks  in  Fig.  (4b)  are  associated  with  a  steep  slope  in  the  integral  curves.  The  set  grouping  is  given  in  Table  I. 
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Radial  profile  with  axial  averaging 


Shot  1484  pinhole  photograph 


r  (cm) 


Shot  1 498  pinhole  photograph 


Figure  6.  Two  illustrations  of  the  method  of  determining  the  width  of  the  radiating  region  by  averaging  the  film  darkness 
in  the  axial  (z)  direction.  Radial  coordinates  given  here  measure  distance  on  the  film,  which  is  twice  the  actual  experimental 
distance.  At  the  top  is  the  pinhole  photograph  and  corresponding  radial  profile  for  shot  1484,  a  fairly  well-defined,  straight 
implosion.  When  averaged  in  z,  the  pixel  values  (black=0,  white=255)  form  a  well-defined,  comparatively  narrow  peak, 
from  which  a  "pinch  radius"  can  be  obtained.  At  the  bottom  is  the  same  information  for  shot  1498,  an  implosion  which  is 
still  fairly  well-defined,  but  which  has  noticable  bends  and  hot-spots.  For  this  shot,  the  pinch  profile  is  broad,  and  has  two 
peaks,  making  determination  of  the  pinch  boundaries  much  more  difficult. 
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Pinhole  photograph  for  shot  1 606 


Figure  7a.  Pinhole  photograph  and  analysis  from  shot  1606,  having  similar  parameters  to  shot  1483  (for  which 
there  was  no  pinhole  data).  Shown,  as  a  function  of  axial  position  z  are  1)  pinhole  photograph;  2)  pinch  diameter, 
determined  as  described  in  the  text;  3)  image  contrast  (maximum  255);  4)  position  of  darkest  point  ("center"  of 
pinch).  The  overall  pinch  width  was  taken  to  be  the  average  of  the  values  in  (2),  except  that  points  were  not 
included  when  contrast  was  less  than  50.  The  center-point  position  gives  some  indication  of  pinch  straightness. 
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Pinhole  photograph  for  shot  1484 
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Figure  7b.  Pinhole  photograph  and  analysis  from  shot  1484.  Shown,  as  a  function  of  axial  position  z  are 
1)  pinhole  photograph;  2)  pinch  diameter,  determined  as  described  in  the  text;  3)  image  contrast  (maximum  255); 
4)  position  of  darkest  point  ("center"  of  pinch).  The  overall  pinch  width  was  taken  to  be  the  average  of  the  values 
in  (2),  except  that  points  were  not  included  when  contrast  was  less  than  50.  The  center-point  position  gives  some 
indication  of  the  straightness  of  the  pinch. 
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Pinhole  photograph  for  shot  1488 


Figure  7c.  Pinhole  photographs  and  analysis  from  shot  1488.  Shown,  as  a  function  of  axial  position  z  are 
1)  pinhole  photograph;  2)  pinch  diameter,  determined  as  described  in  the  text;  3)  image  contrast  (maximum  255); 
4)  position  of  darkest  point  ("center"  of  pinch).  The  overall  pinch  width  was  taken  to  be  the  average  of  the  values 
in  (2),  except  that  points  were  not  included  when  contrast  was  less  than  50.  The  center-point  position  gives  some 
indication  of  the  straightness  of  the  pinch. 
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•Figure  7d.  Pinhole  photograph  and  analysis  from  shot  1489.  Shown,  as  a  function  of  axial  position  z  are 
1)  pinhole  photogr^h;  2)  pinch  diameter,  determined  as  described  in  the  text;  3)  image  contrast  (maximum  255); 
4)  position  of  darkest  point  ("center"  of  pinch).  The  overall  pinch  width  was  taken  to  be  the  average  of  the  values 
‘in  (2),  except  that  points  were  not  included  when  contrast  was  less  than  50.  The  center-point  position  gives  some 
indication  of  the  straightness  of  the  pinch. 
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Figure  7e.  Pinhole  photograph  and  analysis  from  shot  1498.  Shown,  as  a  function  of  axial  position  z  are 
1)  pinhole  photograph;  2)  pinch  diameter,  determined  as  described  in  the  text;  3)  image  contrast  (maximum  255); 
4)  position  of  darkest  point  ("center"  of  pinch).  The  overall  pinch  width  was  taken  to  be  the  average  of  the  values 
in  (2),  except  that  points  were  not  included  when  contrast  was  less  than  50.  The  center-point  position  gives  some 
indication  of  the  straightness  of  the  pinch. 
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Figure  8a.  Pinhole  photogr^h  and  analysis  from  shot  1492.  Shown,  as  a  function  of  axial  position  z  are 
1)  pinhole  photograph;  2)  pinch  diameter,  determined  as  described  in  the  text;  3)  image  contrast  (maximum  255); 
4)  position  of  darkest  point  ("center"  of  pinch).  The  overall  pinch  width  was  taken  to  be  the  average  of  the  values 
In  (2),  except  that  points  were  not  included  when  contrast  was  less  than  50.  The  center-point  position  gives  some 
indication  of  the  straightness  of  the  pinch. 
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8b.  Pinhole  photograph  and  analysis  from  shot  1493.  This  low-contrast  shot  shows  a  discernible  pinch 
radius,  but  suffers  from  noise  around  z=0.6  cm.  Shown,  as  a  function  of  axial  position  z  are  1)  pinhole  photograph; 

2)  pinch  diameter,  determined  as  described  in  the  text;  3)  image  contrast  (maximum  255);  4)  position  of  darkest 
point  ("center"  of  pinch).  The  overall  pinch  width  was  taken  to  be  the  average  of  the  values  in  (2),  except  that  points 
were  excluded  when  contrast  was  less  than  50.  The  center-point  position  gives  some  indication  of  pinch  straightness. 
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'Figure  8c.  Pinhole  photograph  and  analysis  for  shot  1494.  Low  contrast  favors  computer  determination  of  pinch 
radius.  Shown,  as  a  function  of  axial  position  z  are  1)  pinhole  photograph;  2)  pinch  diameter,  determined  as 
described  in  the  text;  3)  image  contrast  (maximum  255);  4)  position  of  darkest  point  ("center"  of  pinch).  The 
‘overall  pinch  width  was  taken  to  be  the  average  of  the  values  in  (2),  except  that  points  were  not  included 
when  contrast  was  less  than  50.  The  center-point  position  gives  some  indication  of  the  straightness  of  the  pinch. 
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Figure  8d.  Pinhole  photograph  and  analysis  for  shot  1495.  Shown,  as  a  function  of  axial  position  z  are 
1)  pinhole  photograph;  2)  pinch  diameter,  determined  as  described  in  the  text;  3)  image  contrast  (maximum  255); 
4)  position  of  darkest  point  ("center"  of  pinch).  The  overall  pinch  width  was  taken  to  be  the  average  of  the  values 
in  (2),  except  that  points  were  not  included  when  contrast  was  less  than  50.  The  center-point  position  gives  some 
indication  of  the  straightness  of  the  pinch. 


30 


Figure  8e.  Pinhole  photograph  and  analysis  for  shot  1496.  This  photograph  was  obviously  ruined  by  hght 
leakage  on  the  top  (r=0).  Shown,  as  a  function  of  axial  position  z  are  1)  pinhole  photogr^h;  2)  pinch  diameter, 
determined  as  described  in  the  text;  3)  image  contrast  (maximum  255);  4)  position  of  darkest  point  ("center"  of 
pinch). 


Figure  8f.  Pinhole  photograph  and  analysis  for  shot  1497.  Shown,  as  a  function  of  axial  position  z  are 
1)  pinhole  photograph;  2)  pinch  diameter,  determined  as  described  in  the  text;  3)  image  contrast  (maximum  255); 
4)  position  of  darkest  point  ("center"  of  pinch).  The  overall  pinch  width  was  taken  to  be  the  average  of  the  values 
in  (2),  except  that  points  were  not  included  when  contrast  was  less  than  50.  The  center-point  position  gives  some 
indication  of  the  straightness  of  the  pinch. 
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Radius  of  radiating  region 
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Figure  9.  Results  of  the  pinhole  photograph  analysis.  The  measured  radius  of  the  radiating  region 
(by  the  FWQM  method),  plotted  as  a  function  of  array  linear  mass  density  and  of  initial  array 
radius.  Measurements  from  set  1  axe  plotted  with  a  plus  (+),  while  those  from  set  2  are  plotted  with 
a  star  (*). 
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Figure  10.  Results  of  CRE  theoretical  calculations  showing  the  values  of  electron  temperature  and 
ion  density  that  would  generate  various  values  of  keV  radiative  power,  for  a  pinch  of  radiating 
diameter  3  mm  (e.g.,  shot  1483). 
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Rgure  11.  Results  of  CRE  theoretical  calculations  showing  the  values  of  electron  temperature 
and  ion  density  that  would  generate  various  values  of  the  Ly  a  to  He  o  line  ratio,  for  a  pinch  of 
radiating  diameter  3  mm  (e.g.,  shot  1483).  For  these  contours,  the  He  a  line  strength  includes  the 
intercombination  line. 
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Figure  13.  Ion  density,  inferred  from  the  spectroscopic,  pulsewidth  and  pinhole  measurements, 
plotted  as  a  function  of  both  array  linear  mass  density  and  of  initial  array  radius.  Measurements 
from  set  1  are  plotted  with  a  plus  (+),  while  those  from  set  2  are  plotted  with  a  star  (*). 
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Figure  14.  Electron  temperature,  inferred  from  the  spectroscopic,  pulsewidth  and  pinhole 
measurements,  plotted  as  a  function  of  both  array  linear  mass  density  and  of  initial  array  radius. 
Measurements  from  set  1  are  plotted  with  a  plus  (+),  while  those  from  set  2  are  plotted  with  a  star 
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Figure  15.  Fraction  of  the  total  imploding  mass  that  is  located  in  the  radiating  region.  Computed 
from  the  measured  pinch  radius,  the  inferred  ion  density,  and  the  known  wire  array  mass.  The  mass 
fraction  is  plotted  as  a  function  of  both  the  initial  linear  mass  density  of  the  array  and  of  initial  array 
radius.  Measurements  from  set  1  are  plotted  with  a  plus  (+),  while  those  from  set  2  are  plotted  with 
a  star  (*). 
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Figure  16.  The  keV  radiative  power,  as  a  function  of  both  the  parameter  rj,  which  indicates  the 
average  kinetic  energy  per  particle  at  stagnation;  and  of  the  fraction  of  total  mass  in  the  radiating 
region  (plotted  in  Fig.  (15)).  Measurements  from  set  1  are  plotted  with  a  plus  (+),  while  those  from 
set  2  are  plotted  with  a  star  (*). 
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